Isospectrality in Chaotic Billiards 
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We consider a modification of isospectral cavities whereby the classical dynamics changes from 
pseudointegrable to chaotic. We construct an example where we can prove that isospectrality is 
retained. We then demonstrate this explicitly in microwave resonators. 



PACS numbers: 05.45.Ac, 03.65. Ge, 41.20.-q 

Recently it has been shown that it is possible to con- 
struct two drums which have different shapes but sound 
exactly the same [1]. This answers the famous question 
asked by Kac [2] in 1966 "Can you hear the shape of 
a drum?" the answer being "no" . Gordon, et. al. con- 
structed an example of a pair of two-dimensional domains 
which had different shapes but had identical eigenvalue 
spectra for the Laplace operator. Since then, a large 
number of such isospectral pairs have been obtained. 

One common feature of all shapes constructed so far is 
that they are mostly polygonal. Hence the classical dy- 
namics of a particle in billiards of these shapes is pseudo- 
integrable. A question of interest then is whether isospec- 
trality can be achieved even for cavities with chaotic dy- 
namics, which is typical of domains that have convex 
pieces, and hence are non-polygonal. We address this 
question, viz. are there sound-alike chaotic drums, both 
theoretically and through experiments using microwave 
resonators. 

Isospectrality is fundamentally a consequence of topol- 
ogy. The essential aspects of isospectrality can be proved 
using the example of the two isospectral domains CI 
and C2 shown in Fig. 1. The proof consists in show- 
ing that given any eigenfunction in one domain we can 
construct a corresponding one in the other domain, with 
the same eigenvalue, and vice-versa. Each domain con- 
sists of seven distinct sub-domains each in the form of 
a triangle. We label these sub-domains in an arbitrary 
fashion, using numbers 1, 2... 7 for domain CI and the al- 
phabets A, B, ...G for domain C2. Note that edges of the 
triangles are marked differently (By dotted, dashed and 
solid lines) and this allows us to make a unique corre- 
spondence between any pairs of triangles. Consider any 
wavefunction, 0, in domain CI, which satisfies the eigen- 
value equation — V 2 = k 2 ip with Dirichlet boundary 
conditions (0 vanishes on the boundary of the domain). 
Let us denote by 0i the restriction of the wavefunction 
in sub-domain i (i.e. ipi(r) — 0(r) if r is a point in 
the ith sub-domain, else ipi{f) = 0). Similarly we can 
define the restricted wavefunctions {0a, V'b---V , g} from 
any wavefunction in domain C2. Starting from the wave- 
function -0 in CI let us construct the following restricted 
wavefunctions in domain C2: 

IpA = 02 ~ 01 + 07 



4>B 
1p c = 

IpD 
IpE 

i>G = 



03 + 01 + 05 
"03 + 02+04 

04 - 01 + 06 

05 - 02 - 06 
07 - 03 + 06 

-07 + 05 -04 
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The notation used requires some explanation: to con- 
struct 0,4 = 02 — 01 + 07, we first move the three 
domains 1, 2 and 7 so that they are on top of each 
other and all similarly marked edges coincide. This 
may require us to flip domains about one of the bases 
and in such cases we have denoted the wavefunction 
with a tilde (e.g. 0i). The wavefunction 0a is then 
obtained by adding (or subtracting) the values of the 
three functions at each point. It is easy to see that 
0' = ipA + 0s + 0c + 0n + 0.E + 0_f + 0g is an eigenfunc- 
tion for domain C2 with the same eigenvalue. For this we 
notice that: 

(1) Laplace's equation is satisfied in every domain and 

(2) It can be verified that the wavefunction vanishes on 
the boundary and matches smoothly across sub-domains. 
For example consider the sub-domains A and B which are 
separated by a dashed line. The wavefunctions arc given 
by 02 — 0i + 07 and ips = 03 + ipi + 05 • The smoothness 
follows since from the wavefunction in CI we see that 02 
matches smoothly with 03 across the dashed boundary, 
similarly 07 matches 05 and —tpi matches 0i. 

Similarly one can construct an eigenfunction for CI 
starting from any given eigenfunction in C2. Thus we 
have demonstrated a one-to-one correspondence between 
the states in the two cavities and hence proved isospec- 
trality. 

We now modify the domain geometry so as to make 
the dynamics chaotic. It is expected that making a part 
of the boundary convex (inwards into the domain) should 
make the dynamics chaotic. This is related to the fact 
that on such boundaries, any two particle trajectories 
which are close to each other, diverge rapidly after being 
reflected [5]. A well known example of a chaotic billiard 
is the Sinai billiard obtained by placing a circular scat- 
terer inside a square. In our case, to obtain the modified 
geometry we first place a scatterer of arbitrary shape in- 
side one of the triangular sub-domains of any one domain 
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and then place one in a similar position in every other tri- 
angle. An example with disc-shaped scatterers is shown 
in Fig. 2. The identification of edges on the two domains 
makes this construction unique. Thus notice that, in ev- 
ery triangle, the scatterer is placed close to a vertex where 
a solid and dotted edge meet. The wavefunction in each 
domain now changes since it has to vanish on and inside 
the boundary of the scatterers. Our construction of the 
modified geometry with scatterers is such that the proof 
for isospectrality given above can be repeated, since the 
wavefunctions still satisfy the relations given by eq.(l). 

A direct physical proof of isospectrality can be ob- 
tained by experiments utilizing microwave cavities [6,7] 
which provide a simple and powerful method of simulat- 
ing single particle time-independent quantum mechanics 
in two dimensions. This follows from the fact that un- 
der appropriate geometrical constraints, Maxwell's equa- 
tions in a cavity reduces to the Schrodinger equation of 
a free particle inside a two dimensional domain of arbi- 
trary shape and topology . Infact one can show that, for a 
cavity with thickness (in the z direction, say) small com- 
pared to the dimensions in the other transverse directions 
(in the xy-plane) , the z-component of the microwave elec- 
tric field ^(x^y) = E z satisfies the time- independent 
Schrodinger wave equation — {d 2 + d 2 )^ = k 2 ^ (with 
the identification k = 2wf/c, f being the frequency and 
c the speed of light) and ip vanishes on the boundary 
of the domain. This correspondence is exact for all fre- 
quencies / < c/2d where d is the thickness of the cavity. 
We note that this is also the Helmholtz equation which 
describes, for example, vibrations of a drum. Using this 
equivalence, various phenomena (such as quantum chaos) 
have been studied. Isospectrality has earlier been demon- 
strated by Sridhar and Kudrolli [8] using microwave cav- 
ities shaped as in Fig. 1. In the experiments one obtains 
the resonance modes of the cavities. Thus the microwave 
transmission spectra directly yields the eigenvalues of the 
cavity being measured. The advantage of this approach 
is that it can be easily applied to arbitrary 2-D domains, 
for which numerical simulations are very hard [9,10] and 
may sometimes be practically impossible. 

In our experiments we consider the same set of cavities 
(Fig. 1) as the ones considered by Sridhar and Kudrolli 
[8] and investigate the question of isospectrality in the 
presence of scatterers placed in the specified way inside 
the cavities. A schematic of the experimental cavity is 
shown in Fig. 3. The desired domain is cut out from 
a brass plate of thickness d = 6mm. Two other brass 
plates are placed on top and below the hole to form a 
closed cavity. As shown in the figure microwaves were 
coupled in and out using loops terminating coaxial lines 
that enter through the sides of the cavity. The length 
of bases of the triangular sub-domains was taken to be 
a = 8cm and the thickness of the cavity was d = 6mm. 
The small thickness of the cavity makes it essentially 2- 
dimcnsional and the correspondence between Maxwell's 
equations and Schrodinger's equation is good for frequen- 
cies / < c/2d — 25GHz. For all metallic objects in the 



2-D space between the plates, Dirichlet boundary condi- 
tions apply inside the metal. 

All measurements were carried out using an HP8510B 
vector network analyzer which measured the transmis- 
sion (S*2i) parameters. The typical values of quality fac- 
tor obtained range from a maximum of 850, at the lower 
end of the spectrum, to a minimum of 250. 

Results: We first attempt to reproduce the results 
in [8] for the cavities shown in Fig. 4. We show in 
Fig 4 the traces of the spectrum for the two cavities 
in the frequency range 1 — 5 GHz. The first 30 res- 
onances of the two cavities are listed in Table I and 
one sees that the eigenvalues match to better than 1%. 
One sees that each resonance present in one is present 
in the other. A few lines are missing and this is at- 
tributed to the fact that the particular coupling posi- 
tions we used may not excite some modes. The re- 
maining inaccuracies are due to imperfections in the ma- 
chining, and in the clamping together of various parts 
of the cavity. Note, that the amplitudes themselves 
may be different, as that depends on the location of the 
coupling and hence to the way the modes are excited. 

TABLE I. The table lists the first 33 resonances in the two 



cavities. 


Resonant frequency 


Resonant frequency 


percentage 


in CI (in MHz) 


in C2(in MHz) 


discrepancy 


1902.500 


1903.750 


0.0657 


2271.250 


2274.375 


0.1374 


2700.625 


2719.375 


0.6895 


3045.625 


3062.500 


0.5510 


3217.500 


3215.000 


-0.0778 


3612.500 


3631.250 


0.5163 




3892.500 




4054.375 






4184.375 


4200.625 


0.3868 


4303.125 


4328.125 


0.5776 


4488.125 


4528.125 


0.8834 


4743.750 


4756.250 


0.2628 


4898.750 






5026.875 


5031.875 


0.0994 


5171.875 


5194.325 


0.4322 


5426.250 


5474.325 


0.8782 


5488.750 






5625.625 


5627.500 


0.0333 


5793.750 


5808.750 


0.2582 




5903.125 




5928.750 


5940.625 


0.1999 


6085 






6222.500 


6242.500 


0.3204 


6312.500 


6352.500 


0.6297 


6497.500 


6492.500 


-0.0770 


6680.000 


6707.500 


0.4100 


6750.000 


6775.000 


0.3690 




6790.000 




6855.000 


6877.500 


0.3272 


6930.000 


6992.500 


0.8938 
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Thus we have obtained the energy spectrum for the given 
set of isospectral cavities and verified that each eigen- 
value in one is present in the other at the same resonance 
value. 

As a check on the quality of the spectral data we 
compare the cumulative number of resonance levels as 
a function of frequency, obtained experimentally, with 
the Weyl formula for the integrated density of states in 
a two-dimensional domain [11]: 

, -r- / % Ak Sk 

n w = -*f-*; +k > (2) 

where A and S are the area and perimeter of the do- 
main, and K is a correction term associated with its 
topology. For a polygonal billiard with inner angles 
a t this is given by K = ~ T")- In thc 

present case we find K = 0.42. We show in Ta- 
ble II a comparison between the experimental results 
with the above formula. The agreement is quite good. 



TABLE II. The table gives a comparison of the cumulative 
frequency in the cavities with the Weyl estimate 



frequency 


cumulative 


Weyl's estimate 


(in GHz) 


resonant frequency 




1 








2 


1 


0.8 


3 


3 


3.4 


4 


7 


7.5 


5 


13 


13.1 


6 


21 


20.4 


7 


30 


29.2 


8 


38 


39.5 



Results for the chaotic geometry: We place the scatter- 
ers inside the cavity following the prescription outlined 
above. The scatterers are taken to be metallic cylinders 
of diameter 1.0 cm and height equal to the thickness of 
the cavities. The modified spectrum from the two cavi- 
ties is shown in Fig. 5 and we list the first 22 resonances 
in Table III. We again find that the eigenvalues in the 
two cavities match to within 1%. Thus there is clear 
evidence that isospectrality is retained in the modified 
chaotic geometry. 

It may be noted that thc introduction of thc scat- 
terers changes the topology of the domain from be- 
ing simply-connected to now being multiply-connected. 
We now have a polygonal box with p = 7 circular 
holes. In this case the topology term in Weyl's formula 
for the integrated density of states is given by: K = 
53(tJ--^)-|- ^ e comparison with Weyl's formula 
is given in Table IV. The agreement is not very good at 
low frequencies and the number of levels seems to be 
somewhat higher than that given by the Weyl estimate. 



TABLE III. The first 22 resonant frequencies in the cavities 



with scatterers 


Resonant frequency Resonant frequency 


percentage 


in CI (in MHz) 


in 


C2(in MHz) 


discrepancy 


2181.875 




2175.000 


-0.3161 


2330.000 




2338.875 


0.3795 


2906.250 




2933.125 


0.9163 


3303.750 




3311.250 


0.2265 


3346.875 




3361.875 


0.4462 


3746.875 




3766.250 


0.5144 


4135.625 




— - 


— - 


4175.625 




4175.625 


0.000 


4207.500 




4222.500 


0.3552 


4645.000 




4666.250 


0.4554 


4966.875 




4956.250 


-0.2144 


5056.250 




— 


— - 


5106.250 




5101.250 


-0.0980 


5148.750 




5124.375 


-0.4757 


5540.000 




5576.875 


0.6612 


5603.125 




5633.125 


0.5326 






5889.375 




6045.000 




6002.500 


-0.7080 


6435.000 




6450.000 


0.2325 


6490.000 




6505.000 


0.2306 


6850.000 




6872.500 


0.3274 


6955.000 




6995.000 


0.5718 


TABLE IV. 


Comparison of Weyl's estimate and cumula- 


tive frequency obtained experimentally after placing the scat- 


terers 








frequency 


cumulative 


Weyl's estimate 


(in GHz) 


resonant frequency 




1 










2 










3 


3 




0.9 


4 


6 




4.5 


5 


11 




9.7 


6 


17 




16.3 


7 


22 




24.5 



To make the demonstration more convincing and illus- 
trate the non-triviality of the isospectral construction 
with scatterers, we consider another geometry (Fig. 6) 
where the scatterers in the second cavity are placed in a 
somewhat different manner. The arrangement still seems 
to follow the folding construction and naively one would 
expect isospectrality. However on closer inspection one 
finds that the correct correspondence between the edges 
of the sub-domains has not been satisfied and the wave- 
function matching condition infact no longer holds and 
so we should not get isospectrality. We plot the spectrum 
for this case in Fig. 7. We see a marked difference from 
that in Fig. 5, namely we find that there is no correspon- 
dence between the spectral lines from the two cavities. 
This shows clearly that isospectrality is indeed obtained 
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only for the special arrangement of scatterers in Fig. 2. 

In conclusion wc have demonstrated that isospectral- 
ity is unrelated to the underlying classical dynamics of 
a particle. We have shown a simple way of introducing 
scatterers of arbitrary shape into polygonal cavities in 
such a way that isospectrality is retained. This leads us 
to a new class of isospectral scatterers and also a bet- 
ter understanding of the essential features necessary for 
isospectrality. 
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FIG. 1. Isospectral cavities CI and C2. The outer edges of 
the polygonal structure constitute the boundary of the cav- 
ity. The inner edges have been marked to show the seven 
triangular sub-domains within each cavity. 




FIG. 2. Isospectral cavities with scatterers in the shape of 
discs. The wavefunction vanishes on the boundary and inside 
of every scatterer. 




FIG. 3. A schematic of the experimental cavity. This shows 
a brass plate on which a hole of the desired cavity shape has 
been cut. This plate is sandwiched between two other brass 
plates to form a closed cavity. 



4 





•5 



